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Abstract 

Using linear perturbation theory and the Priedmann-Lemaitre solutions of the cosmological 
field equations, we derive analytically a second-order differential equation for the evolution of the 
linear bias factor, b{z), between the background matter and a mass-tracer fluctuation field. We 
find b{z) to be a strongly dependent function of redshift in all cosmological models. Comparing 
our analytical solution with the semi-analytic model of Mo & White, which utilises the Press- 
Schechter formalism and the gravitationally induced evolution of clustering, we find an extremely 
good agreement even at large redshifts, once we normalize to the same bias value at two different 
epochs, one of which is the present. Furthermore, our analytic b{z) function agrees well with the 
outcome of N-body simulations. 
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1. Introduction 

The concept of biasing between different classes of extragalactic objects and the background 
matter distribution was put forward by Kaiser (1984) and Baxdeen et al. (1986) in order to explain 
the higher amplitude of the 2-point correlation function of clusters of galaxies with respect to that 
of galaxies themselves. 

In this framework biasing is assumed to be statistical in nature; galaxies and clusters are 
identified as high peaks of an underlying initially Gaussian random density field. Biasing of galaxies 
with respect to the dark matter distribution was also found to be an essential ingredient of CDM 
models of galaxy formation in order to reproduce the observed galaxy distribution (cf. Davies et 
al. 1985; Benson et al. 2000). 

The classical approach to study the rcdshift evolution of bias utilises the ratio of the correlation 
functions of objects and dark matter, which are assumed to be related via the square of a scale 
independent bias factor. However, in this study we will use the definition by which the extragalactic 
mass tracer (galaxies, halos, clusters) fluctuation field, Str, is related to that of the underlying mass, 
Sm, by 

Str = bSra , (1) 

where b is the linear bias factor. Note that the former definition results from the latter but the 
opposite is not necessarily true. The bias factor may have many dependencies; even assuming that 
it is scale independent, it necessarily depends on the type of the mass tracer as well as on the epoch 
z, since the fluctuations evolve with time as gravity draws together galaxies and mass. It is evident, 
therefore, that the bias factor should also depend on the different cosmological models and dark 
matter content of the Universe (for a recent overview see Klypin 2000). 

More realistic biasing schemes have been proposed in the literature. Coles (1993) introduced 
the idea of biased galaxy formation in which galaxies form with a probability given by an arbitrary 
function of the local mass density. Mann, Peacock & Heavens (1998) investigated the properties of 
different bias models of galaxy distributions that results from local transformations of the present- 
day density field. The deterministic and linear nature of eq.(l) has been challenged (cf. Dckel & 
Lahav 1999; Tegmark &; Bromley 1999) and indeed some non-linearity of the biasing relation is 
necessary to reconcile high biasing with deep voids. Despite the above, the linear biasing assumption 
is still a useful first order approximation which, due to its simplicity, it is used in most studies of 
large scale (linear) dynamics (cf. Strauss 8z Willick 1995 and references therein; Branchini et al. 
1999; Schmoldt et al. 1999; Plionis et al. 2000). 

Different studies have indeed shown that the bias factor is a monotonically decreasing function 

of redshift. An important advancement in the analytical treatment of the bias evolution was the 
work of Mo &: White (1996) in which they used the Press- Schcchter (1974) formalism and found 
that in an Einstein-de Sitter universe the linear bias factor evolves strongly with redshift. Using a 
similar formalism, Matarrese et al. (1997) extended the Mo &; White results to include the effects 
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of different mass scales (see also Catelan et al 1998). 

Steidel et al. (1998) confirmed that the Lyman-break galaxies are very strongly biased tracers 
of mass and they found that b{z = 3.4);^ 6, 4, 2, for SCDM, ACDM (O = 0.3) and OCDM (O = 0.2), 
respectively (sec also Giavalisko et al 1998). A similar value for the ACDM model was obtained by 
Cen &: Ostrikcr (2000) using high resolution Nbody/hydro simulations in which they treated DM, 
gas as well as star formation. The use of high resolution N-body simulations (cf. Klypin et al. 1996; 
1999, Cole et al. 1997 and references therein) have shown that anti-biasing {b < 1) should exist at 
scales r ~ 3 — 8h~^ Mpc, for the open and flat low-0 models, in contrast with = 1 models, where 
b > 1. Cohn et al (1999), using high-resolution N-body simulations of SCDM, ACDM, OCDM 
and rCDM models, which avoid the so called "overmerging" problem, found that indeed biasing 
evolves rapidly with redshift, while Kauffmann et al. (1999) combining semi-analytic models of 
galaxy formation and N-body simulations has also studied the evolution of clustering in different 
cosmologies. 

In this paper we will not indulge in such aspects of the problem but rather, working within 
the paradigm of linear and scale-independent bias, we will derive the functional form of its redshift 
evolution in the matter dominated epoch and in all cosmological models. The Einstein de-Sitter 
case has been studied in the past (cf. Nusser &; Davis 1994; Pry 1996; Bagla 1998) using the 

continuity equation, which is a first order differential equation, to derive a solution, oc (l + z), valid 
only for low z's. Our approach is to use the perturbation evolution equation which combines the 
continuity, the Euler and the Poisson equations and which is a second order differential equation. 
We should therefore expect to find a further component to the known solution. 

The paper is organised as follows: in section 2 we discuss the basic models for the linear bias 
evolution, in section 3 we derive the basic differential equation describing the evolution of the linear 
bias factor, while in section 4 we present its analytical solution for the different cosmological models 
and a comparison with previous models and N-body simulation results. Finally, in section 5 we 
summarise our main results. 

2. Models for bias evolution 

Theoretical expectations regarding the cosmological evolution of bias have been investigated 
using analytical calculations, semi-analytical approximations and N-body simulations. In this sec- 
tion we shortly describe some of these models in order to compare them with our results. 

2.1. Test Particle or Galaxy Conserving Bias (Ml): 

This model, proposed by Nusser & Davis (1994), Fry (1996), Tegmark & Peebles (1998), 
predicts the evolution of bias, independent of the mass and the origin of halos, assuming only that 
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the test particles fluctuation field is related proportionally to that of the underlying mass. Thus, 
the bias factor as a function of redshift can be written: 

b{z) = 1 + (5o - l)D{z)-^ (2) 

where bo is the bias factor at the present time. Bagla (1998) found that for SCDM model and in 
the range < 2; < 1 the above formula describes well the evolution of bias. 



2.2. Halo Models (M2): 



Mo &: White (1996) using the Press-Schecter formalism, have developed a model for the evo- 
lution of the correlation bias, which depends on halo mass, and found, in an Einstein-de Sitter 
Universe, that: 



biz) = 1 + 



1 



D{z) 



(3) 



with a reference redshift ^ , Sc the critical overdensity for a spherical top-hat collapse model and 
D{z) = (1 + z)~^ is the linear growth rate of clustering. Parametrising this equation to the present 
epoch one gets: 

b{z) = 0.41 + (6(0) - OAl)D{z)-^ (4) 

Similarly, Matarrcsc ct al. (1997) parametrising the evolution of bias for halos above a certain mass 
M, obtain a similar expression for an Einstein-de Sitter Universe: 



b{z) = 0.41 + (6eff(0) - OAl)D{z)-^ 



(5) 



with 5efr the bias of a sample of halos with a range of masses and /3 depending on the minimum 
mass scale that contributes to the halo correlation function (with (3 <2). 



3. Basic Equations 

The central issue here is to derive the basic differential equation which describes the evolution 
of bias. The present analysis is based on linear perturbation theory in the matter dominated epoch 
(cf. Peebles 1993) and it is an extension of the Ml model. 



^Consider the distribution of halos of mass M, or larger, before typical halos of this mass have collapsed. Therefore, 
to quantify this, the parameter i'{M,z) = 5c/cj{M,z), which related directly with the bias, 6=1 + — 1), is 

usually utilized (where 5c = 1-69 is the critical overdensity for spherical collapse at z = 0, and cr(M, z) is the rms 
linear mass fluctuation on the scale M of halos linearly extrapolated to redshift z) . Thus, z» is fixed by requiring 
v{M,z,) = 1 (cf. Bagla 1998). 
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The time evolution equation for the mass density contrast, (5m = {5p/p)ra, modelled as a 
pressureless fluid with general solution of the growing mode: 6^ ~ A{x)D{t), is (cf. Padmanabhan 
1993): 

+2Hit)^=A7:Gp^S^ , (6) 



dt 



Assuming for simplicity that the mass tracer population is conserved in time, ie., that the effects 
of non-linear gravity and hydrodynamics (merging, feedback mechanisms etc) do not significantly 
alter the population mean, then a similar evolution equation, containing in the right hand side 
the gravitational contributions of all the perturbed matter, should be satisfied for 6tr (see also Fry 
1996; Catelan et al. 1998): 



' dt 

Differentiating twice eq.(l) and using eq.(7) and eq.(6) we obtain: 



bdra + 2 



b + 



Then from equation eq.(8), eq.(6) and Sm — A{x)D{t) we have: 

bD{t) + 2 lolt) + H{t)D{t) \ b + ATTGpraD{t)b = AnGpn,D{t) . 

In order to transform eq.(9) from time to redshift we use the following expression: 

1 



dt 
dz 

where the Hubble parameter is given by: 



with 



H{z) = HoE{z) 



E{z) = [^^(l + zf + ^nil + z)2 + J^a] 



1/2 



(7) 



(8) 



(9) 



(10) 



(11) 



(12) 



and Q, = S'KGpo/'iHg (density parameter), ^Ir = {HgaoR)^^ (curvature parameter), Oa = A/3i?2 
(cosmological constant parameter) at the present time which satisfy Q, + Clji + J^a = 1 and Ho is 
the Hubble constant. 



Finally, the growing solution (cf. Peebles 1993) as a function of redshift is: 



D{z) 



5nE{z) 



jl + x) 



dx 



(13) 



Therefore, as the time evolves with redshift, utilised eq.(lO), eq.(ll) eq.(12) and the relation 



47rGpni = 4:TrGpo{l + z) 



3-^^l + .)3, 



(14) 
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then the basic differential equation for the evolution of the linear bias parameter takes the following 
form: 

d'^b db 

-P{z)- + Q{z)b = Q{z) (15) 



with basic factors, 



and 



dz^ dz 

P(z) = ^ ^d^_^d^ 

^ ^ 1 + z E{z) dz D{z) dz ^ ' 

^ m{l + z) 

' 2EHz) ■ ^ ' 



It is obvious that the above generic form depends on the choice of the background cosmology. Thus, 
the functional form which satisfies the general bias solution for all of the cosmological models is: 

b{z) = y {z; n;n a) + 1 (18) 

where y is the general solution of the homogeneous differential equation: 

g_P(,)^+Q(,)y = . (19) 

Whereas it is obvious that the present theoretical approach takes into account the gravity field, it 
does not interact directly with the nature of the DM particles. 



4. Bias Evolution in different Cosmological Models 

In this section using both eq.(19) and Friedmann-Lemaitre solutions of the cosmological field 
equations we present the analytical solution of bias evolution for the Einstein-de Sitter, the open 
and low-density flat cosmological models. 



4.1. Elements of the Differential Equation Theory 

Without wanting to appear too pedagogical, we remind the reader some basic elements of 
differential equation theory (cf. Bronson 1973). If one is able to flnd any solution yi of eq.(19), 
then a second linearly independent solution be found very easily. Let the second solution can be 
written as 

y2{z) = yi{z)u{z) (20) 

where u{z) is to be determined. Inserting eq.(20) into eq.(19) and remembering that yi satisfies 
the same equation, we find the following equation for u{z): 



dz^ 



+ 



2{dyi/dz) - P{z)yiiz) 

yi 



dz 



(21) 
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Integrating the above equation we have 

dn const 



dz yl{z) 



20 



exp / P{x)Ax (22) 



where zq is an arbitrary initial point. A further integration of eq.(22) yields u{z)^ and inserting 
this value into eq.(20), we obtain the second solution 

y2{z) = const yi{z) [ ^^exp ( / P{t)dt] . (23) 
Jzo Viix) \Jzo ) 

The Wronskian of the two solution y\ and y^ is 

W{yuyi) = yl^ = const exp P{x)dx^ . (24) 

Thus the Wronskian never vanishes which implies that any general solution of eq.(19) is a linear 
combination {y = ciyi + 022/2) of the fundamental set of solutions yi and y2- 



4.2. Einstein - de Sitter Model 

In this case the basic cosmological equations are the following: 

E{z) = il + zf\ (25) 
while the growth factor of the linear density contrast is 

D(^) = (26) 

and thus 

P(z) = (27) 



and 

2{l + z) 



Qi^) = on , .^2 ■ (28) 



It can be found that the function yi{z) = {1 + z) = D''^{z) is a solution of the eq.(19) which is 
to be expected from the Ml model. Therefore, we are looking for the second independent solution 
of the eq.(19). Thus according to the procedure, described before we can calculate the second 
solution directly from eq.(20) y2{z) = (1 + z)^/^ = D~^/'^{z). The general solution of the second 
order differential eq.(19) is the following: 

y{z; 1; 0) = AD-\z) + BD-^'^{z) = A{1 + z) + 8(1 + zf/^ (29) 



with general bias solution b{z) = y{z; 1;0) + 1. To this end this analysis generalise the Ml model 
in the sense that the added function y2{z) dominates the functional form of the bias evolution. Of 
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coursc in order to obtain partial solutions for b{z) we need to estimate the values of the constants 
A and B, which means that we need to calibrate the b{z) relation using two different epochs: 
b{z = 0) = bo and b{z = zi) = bi. Therefore, utilised both the above general bias solution and the 
latter parameters, we can give the expressions for the above constants as a function of bo and bi : 

_ {bo-l)D-'/Hz,)-{b,-l) 

D-3/2(^i)-L>-l(zi) ' ^^^^ 

(6i-l)-(5o-l)iJ-i(^i) 

D-^/^z^) - D-^izi) ■ ^ ' 

For B = (Ml model) we obtain, as we should, ^ = 6o — 1- 

Our generalised solution does not suffer from limitations in the value of b (as does the Ml 
solution); b can take values > 1 and < 1. It is interesting to compare our generalised test-particle 
bias with the more elaborate halo and merging models. Since our approach gives a family of bias 
curves, due to the fact that it has two unknown parameters, (the integration constants A,B), we 
evaluate the latter by using Steidel et al. (1998) value of the bias for Lyman break galaxies which 
gives for ft = 1, 6(3.4) ~ 7. Inserting this into Mo & White (1996) model we obtain 6(0) = 0.75. 
In figure 1 we compare our solution with the Mo & White model and to our surprise we find an 
excellent agreement. This implies that the complete test particle bias solution is an extremely 
good approximation to the more elaborate halo solutions which takes into account, via the Press- 
Schechter formalism, the collapse of different mass halos at the different epochs. 

We further compare our analytic solution with N-body estimates provided by Colin et al. 
(1999) and Kauffmann et al. (1999). In figure 2 our model is represented by lines while the 
numerical results with the different symbols. It is evident that our analytic function, normalized 
to two different epochs of the numerical results, fits extremely well the behaviour of the N-body 
derived bias evolution. 



4.3. Low Density Universes 

The basic cosmological equations in low-density Universes become more complicated than in 
an O = 1 Universe and eq.(19) does not have simple analytical solutions. We therefore present 
approximate analytical solutions which arc valid in the high-rcdshift regime. In order to do so we 
consider that (i) for a low density open Universe, the Einstcin-de Sitter growing mode is a good 
approximation for z'^ — 1 and (ii) for a low density flat Universe the growing mode is well 
approximated by the Einstein de- Sitter case for z^ _ 1 (cf. Peebles 1984b; Carrol, Press & 

Turner 1992). 
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4-3.1. Analytical Approximation for the Open Universe 



Using eqs (16), (17) and z'^ Q, ^ — 1 we obtain the following basic factors of the differential 
equation (19): 

9 O 

(32) 



and 



Piz)- 
Q{z) 



1 + z 2{i + nz) 



(33) 



2{i + z){i + nz) 

where we have used E{z) = (1 + z){l + J7z)^/^ and eq.(26). In this case it can be easily found that 
the function yi{z) = (l + z) + 4(l — 0)/3r2 is a solution of the eq.(19). Thus, after some calculations 
we can obtain the general bias solution: 



b{z) - 1 = y{z;^ (n-^ - iy,n;o) = A 

with 

u{z) = 



{l + z)+4 



(l-») 



+ c 



(1 + 2) + 4 



(1-^) 
3J^ 



{l + zf dz 



(l + z)+4^i^ 



3n 



(1 + i^z)V2 



u{z) (34) 



(35) 



Performing the latter integration one finds that the bias evolution is given by: 



h{z) -1=A 



{l + z)+A 



+ 



2C 



(1 + 2) +4 



(l-») 



(1 - fi) + 3(1 + nz) 

(36) 

It is obvious that for J7 — > 1 the above solution tends to the Einstein-de Sitter case, as it should. 



4.3.2. Analytical Approximations for the A Universe 

In a flat universe with non zero cosmological constant the growing mode approximation leads 
to the following basic factors: 

Piz) = ^ - ^^(^ + ")' (37) 
1 + z 2[J^(1 + 2)3 + Oa] ^^'^ 

and 

where wc have used E{z) = [$7(1 + z)^ + r^A]^''^ and cq.(26). It is obvious that yi{z) = E{z) is a 
solution of eq.(19) and following a similar procedure to that of the previous subsection we obtain 
the general solution: 

h{z) - 1 = y{z-k. {^-^'^ - ^k) = A p{l + zf + + C + zf + ^'^ u{z) (39) 
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where 



{1 + zfdz 



(40) 



The integral of equation (40) is elliptic and therefore its solution, in the redshift range [z, +oo 
can be expressed as a hyper-geometric function. We finally obtain: 



b{z)-i = A[n{i+zf+nA]'/^+^[n{i+zf+nA]^/^{i+z)-^/^F 
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6'2'6' n(l + z) 



■ (41) 



If 0, — > 1 and — ^ the above bias solution tends to the Einstcin-dc Sitter case, as it should. 
Note that for = 0.7 our solution is valid even at low redshifts, since the growing mode of the 
fluctuations evolution is well approximated by the Einstein de Sitter solution. Therefore we present 
in figure 3 a comparison between the results of the high-resolution N-body simulations of Colin et 
al. (1999) and our solution, parametrised to 6(0) = 0.75 and the bias value of Colin et al. (1999) 
at z = 3. As it is evident the agreement is excellent. 



4-3.3. Analytic Solution for two Limiting Cases 

Although, due to the complex form of the P{z) and Q{z) functions in the case of J7 < 1 
Universes, we cannot solve the bias evolution problem analytically for all redshifts, we can produce, 
however, a complete analytical bias evolution solution for the limiting case of O = (Milne Universe) 
or = and n\ = 1 (de Sitter Universe). 

In the former case we have from eqs (16) and (17) that Q{z) = and P{z) = 0. Applying 
these to eq.(19), we find the general bias solution for this special open cosmological model to be: 

b{z) =A{l + z) + B (42) 

Interestingly, also the density fluctuations, S, have such a z-dependence in the = Universe. 

In the de Sitter Universe, which is dominated only by vacuum energy, putting ~ into 
eq.(16) and eq.(17) we obtain again Q{z) = 0, while the P{z) factor is given by: 

P(^) = ^ , (43) 

with solution 

b{z) =A {1 + z)dz + B . (44) 
Performing the latter integration one finds: 

h{z) = Ail + zf + B , (45) 

with initial condition ho = A + B. 
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We know that the two special solutions (42) and (45) do not correspond to realistic Universes. 
Nevertheless, these solutions can operate as limiting cases of the generic problem. For example, 
using the general solution of the bias (eq.l8) and putting Qo = 0, or = 1 and Qo = 0, the 
factors (1 + z) and (1 + z)'^ should survive in the general well. In Figure 4, we compare 

the bias evolution of these special low-density models with the Einstein-de Sitter one, normalising 
them to the same 6(0) and using the results of Steidel et al. (1998). Of course, by no means do we 
imply that these models predict the same value of 6(0) , we only parameterise our solution in order 
to compare the behaviour of b{z) in these limiting cases. 



5. Summary 

We have introduced analytical arguments and approximations based on linear perturbation 
theory and a linear, scale-independent bias between a mass tracer and its underlying matter fluctu- 
ation field in order to investigate the cosmological evolution of such a bias. We derive a second order 
differential equation, the solution of which provides the functional form of the of bias evolution in 
any Universe. For the case of an Einstein-de Sitter Universe, we find an exact solution which is a 
linear combination of the known solution oc (1 + z) (cf. Bagla 1998 and references therein), derived 
from the continuity equation, and a second term oc (1 -|- z)^/^ which dominates. This solution once 
parametrised at two different epochs, compares extremely well with the more sophisticated halo 
models (cf. Mo & White 1996) and with N-body simulations. 

For the two low-density cosmological models wc find exact solutions, albeit only in the high- 
redshift approximation (where the growing mode of perturbations can be approximated by the 
Einstein-de Sitter solution). We also derive analytical solutions for two limiting low-density Uni- 
verses (ie., Q = 0, Oa = and O = 0, Ha = 1). 

We thank Stefano Borgani and Peter Coles for their useful comments and suggestions. We 
also thank the anonymous referee for the critical comments and useful suggestions. M.Plionis 
acknowledges the hospitality of the Astrophysics Group of the Imperial College, were this work was 
completed. 
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Fig. 1. — Comparison of bias evolution between our solution (continuous line) and the functional 
form (eq.4) of the Mo & White (1996) halo model (dots). 
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Fig. 2. — Comparison of bias evolution between our solution and $7 = 1 N-body simulation results, 
(a) Lines represent our model, parametrised to two epochs of the Colin et al. (1999) rCDM and 
SCDM models, while points represent their results, (b) Similarly for the Kauffmann et al. (1999) 
model. 
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Fig. 3. — Comparison of our ACDM {Qa = 0.7) solution (broken line) with the results of high- 
resolution N-body simulations of Colin et al. 1999 (points). 
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Fig. 4. — Evolution of the bias factor for different models parametrised using the Steidel et al. 
(1998) bias values of Lyman-break galaxies at z = 3.4 (6 ~ 6, 4, 2, for SCDM, ACDM (O = 0.3) 
and OCDM (0 = 0.2), respectively). 



